Evolution of the symmetry energy of hot neutron-rich matter formed in heavy-ion 

reactions 
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It is shown that the experimentahy observed decrease of the nuclear symmetry energy with the 
increasing centrality or the excitation energy in isotopic scaling analyses of heavy-ion reactions 
can be well understood analytically within a degenerate Fermi gas model. The evolution of the 
symmetry energy is found to be mainly due to the variation in the freeze-out density rather than 
temperature. The isoscaling analyses are useful for probing the interaction part of the nuclear 
symmetry energy, provided that both the freeze-out temperature and density of the fragments can 
be inferred simultaneously from the experiments. 
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INTRODUCTION 



Information about the symmetry energy of hot 
neutron-rich matter is important for understanding the 
dynamical evolution of massive stars and the supernova 
explosion mechanisms, while the symmetry energy at 
zero temperature is important for determining proper- 
ties of neutron stars at /3-equilibrium. In particular, the 
electron capture rate on nuclei and/or free protons in pre- 
supernova explosions is especially sensitive to the symme- 
try energy at finite temperatures. The electron captures 
drive the collapsing core towards more neutron-rich mat- 
ter. They affect not only the electron degenerate pres- 
sure working against the gravity but also the neutrino 
fluxes carrying away energy from the core0, |M 1^ 
The larger the symmetry energy, the more difficult it is 
for the electron captures to happen. Heavy-ion reactions 
are a unique means to produce in terrestrial laborato- 
ries the hot neutron-rich matter similar to those existing 
in many astrophysical situations. The possibility of ex- 
tracting useful information about the symmetry energy 
from heavy-ion reactions has stimulated much interest in 
the nuclear physics community 0, 0, • Especially, re- 
cent analyses of the isospin diffusion data in heavy-ion 
reactions M lol Il0l| and the size of neutron-skin in ^"^Pb 
in III Ell have already put a stringent constraint on 
the symmetry energy of cold neutron-rich matter at sub- 
normal densities. This has led to a significantly more 
refined constraint on the mass-radius correlation of neu- 
tron stars |l4l including the fastest pulsar discovered very 
recently [l5|. On the other hand, the temperature depen- 
dence of the symmetry energy for hot neutron-rich matter 
has received so far little theoretical attention. 

Among the phenomena/observables identified as po- 
tentially useful probes of the nuclear symmetry energy, 
the isoscaling coefficients of fragments from heavy-ion re- 
actions have been most extensively studied, see e.g., 
E3| for a recent review. Very interestingly, it was found 



recently that the extracted symmetry energy from the 
isoscaling analyses decreases significantly from its stan- 
dard value of about 25 MeV at normal nuclear matter 
density po to much smaller values with the increasing 
excitation energy or centrality in heavy-ion reactions at 
both Fermi |l8l [ijj and relativistic energies mil 113. 
Moreover, an increasing temperature of the fragmenting 
system was found to accompany the decreasing symme- 
try energy in these reactions. However, the fundamental 
origin of this apparent evolution of the symmetry energy 
is still not clear and it is particularly important and in- 
teresting to understand to what degree the evolution is 
due to the density and/or the temperature dependence 
of the symmetry energy. 

In this work, it is shown that the experimentally ob- 
served evolution of the symmetry energy can be well un- 
derstood within a degenerate Fermi gas model at finite 
temperatures. Furthermore, it is found that the evolution 
of the symmetry energy is mainly due to the variation in 
the freeze-out density rather than temperature when the 
fragments are emitted in the reactions carried out under 
different conditions. 



II. NUCLEAR SYMMETRY ENERGY AT 
FINITE TEMPERATURE 



The Equation of State (EOS) of hot neutron-rich mat- 
ter at a temperature T and an isospin asymmetry S = 
{Pn — Pp)/{Pn + Pp) can be written as [23Ll2^ 

E{p, T, 5) = E{p, T,6^0) + E,ym{p, T)S^ + 0(6^). (1) 

The temperature and density dependent symmetry en- 
ergy Esyrnip, T) for hot neutron-rich matter can thus be 
extracted from E^yrr^ip, T) « E{p, T, 5 = 1) - E{p, T,S = 
0). The symmetry energy Esym{p,T) is the energy cost 
to convert all protons in symmetry matter to neutrons at 
the fixed temperature T and density p. For finite nuclei 
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at temperatures below about 3 MeV, the shell structure 
and pairing as well as vibrations of nuclear surfaces are 
important and the symmetry energy was predicted to in- 
crease slightly with the increasing temperature [25l l26l| . 
Interestingly enough, an increase by only about 8% in 
the symmetry energy in the range of T from to 1 
MeV was found to affect appreciably the physics of stel- 
lar collapse, especially the neutralization processes |25j. 
At higher temperatures, one expects the symmetry en- 
ergy to decrease as the Pauli blocking becomes less im- 
portant when the nucleon Fermi surfaces become more 
diffused at increasingly higher temperatures [23. . In 
this work, we use the thermal model of Mekjian, Lee and 
Zamick (MLZ) |23|. While all of our results have also 
been concurrently verified numerically by using the fi- 
nite temperature Hartree-Fock (HF) approach using both 
Skyrme and Gogny forces [23, here we utilize the MLZ 
approach because of its analytically transparent prop- 
erties. The results obtained within the MLZ approach 
are sufficient for the purposes of this work. Our stud- 
ies based on the finite temperature HF calculations will 
be reported elsewhere [2^. The degenerate Fermi gas 
limit of the MLZ thermal model is appropriate for us to 
understand quantitatively the experimentally observed 
evolution of the nuclear symmetry energy. The symme- 
try energy Esym {p, T) has a kinetic contribution and an 
interaction part. In the MLZ model with Skyrme inter- 
actions, the interaction part is temperature independent, 
i.e., 

Esyrn if, T) = E^;Z {p, T) + El^l (p) . (2) 

At low temperatures, it is known that all mean field quan- 
tities are essentially temperature independent 25]. With 
even momentum-dependent Gogny forces, at tempera- 
tures relevant for fragment formation in heavy-ion reac- 
tions, our HF calculations indicate that the interaction 
part Elylj^ is only slightly T-dcpendent. For tempera- 
tures much less than the Fermi energy, T <C Epermi ~ 
36(/9/po)^/^ the kinetic energy per nucleon of a near- 
degenerate two-component Fermi gas is 27] 

E,,^ ip, T, S) = 21.^/3(1 P') + ^^^,{l~ls% (3) 

where u = p/ po is the reduced density. The kinetic part 
of the symmetry energy is thus 

It is interesting to note that the -E^ *"„ (P: ^) decreases 
with —T^ with a rate depending on the density. On the 
other hand, for T ^ Epermi the system becomes a non- 
degenerate classical gas with a kinetic symmetry energy 

of 13 

E'^lriP^T) = ^T[0.177^p- 0.0033(^p)2], (5) 



where A = ^27rh'^/mT is the thermal wavelength of 
nucleons with an average mass m. It is seen that the 
symmetry energy decreases approximately according to 
Esym{p,T) (X p/VT for X'^p ^ 1. The genuine feature 
of a decreasing symmetry energy associated with an in- 
creasing temperature at both the degenerate and non- 
degenerate Fermi gas limits is consistent with predictions 
of the micro scop ic and/or phenomenological many-body 
theories [^[^[^[^[s^. 

At the typical freeze-out temperatures and densities for 
the emission of intermediate mass fragments in heavy- 
ion reactions, the degenerate Fermi gas limit is appro- 
priate except at extremely low densities where cluster- 
ing effects may become important (2£, ^1] . To evaluate 
the relative importance of the T-dependent part AE'^ = 
7r^r^/(1260u^/'^) with respect to the total symmetry en- 
ergy for cold neutron-rich matter, we show in Fig. 1 the 
ratio AE'^ / Esym{p,T = 0). The analyses of the isospin 
diffusion data from NSCL/MSU [ll lj l ^ nd the study 
on the size of neutron-skin in •^°^Pb [ill Il2l ITsI have re- 
cently consistently constrained the symmetry energy of 
cold matter to be around 32(p/po)°'^ < Esym{p,T = 
0) < 32(p/po)^'^ at sub-normal densities. Using the 
above two limits for the Esym {p, T = 0) at a typical 
freeze-out temperature of T = 5 MeV for the fragment 
emission in heavy- ion reactions, it is seen that the ratio 
AE'^ /Esym{p,T = 0) increases quickly with decreasing 
density. The T-dependent part of the symmetry energy 
becomes increasingly more appreciable, e.g., up to about 
35% at p = 0.1/90 for Esym{p,T = 0) = 32(p/po)' ^ 
Moreover, since the ratio increases quadratically with T, 
the effect will be much larger at higher temperatures. 
This result also shows the magnitude, especially at low 
densities and/or high temperatures, of some artificial ef- 
fects that would be introduced should one attribute the 
entire evolution of the symmetry energy to its density 
dependence while neglecting its intrinsic temperature de- 
pendence. 
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FIG. 1: (Color online) The relative importance of the 
temperature-dependent part of the nuclear symmetry energy. 

Corresponding to the above two limits for the 
Esym{p,T = 0), the interaction part of the symme- 
try energy -EJ™(/9) is constrained between E^J^l^{p) = 
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130m - 111.4ui i (labeled as x = 0) and E^l^{p) = 
3.7u + 14.9ui-6 (labeled as x = -1) H |l3- With 
this E"yl^{p) and Eq. (4), we can now examine the 
evolution of the total symmetry energy (Eq. (2)) as 
functions of density and temperature. Shown in Fig. 
2 is the evolution of the normalized symmetry energy 
[Esym{p,T) - Esy„i{po,0)]/ Esy„i{pQ,0) as a function of 



the symmetry energy of the fragments at the freeze-out. 
Moreover, they use the finite size effects, such as the sur- 
face symmetry energy and its temperature dependence, 
to explain the observed smaller value of Csym compared 
to the symmetry energy of about 30 MeV for cold nu- 
clear matter at the saturation density. However, from the 
very nature of the isoscaling phenomenon itself that iso- 
topes/isotones having very different masses (sizes) fall on 
the same curve described by a single scaling coefficient, 
it is hard to believe that the finite-size effects have any 
influence on the Csym at all. In another word, unless the 
finite size effects on both the Csym and the temperature 
T, if there is any at all, are completely cancelled out, the 
isoscaling phenomenon should not have been observed 
in multifragmentation in the first place. Indeed, in the 
AMD analyses of isoscaling in multifragmentation, it was 
found that "the extracted symmetry energy shows almost 
no surface effect in it, which suggests that the properties 
of infinite nuclear matter can be directly obtained from 
the information of fragmentation" j38l | . 

Within the grand canonical model for multifragmen- 
tation, another possible reason for the extracted small 
value of Csym is that fragments themselves at freeze-out 
are dilute. This picture, however, seems to contradict 
the basic Fisher hypothesis that correlations inside a di- 
lute medium are exhausted by clusterization. One possi- 
ble explanation put forward was that primary fragments 
formed in heavy-ion reactions are hot and thus also ex- 
panded [islls^ . However, this explanation is insufficient 
to explain the much smaller value of Csym observed in 
central collisions. Moreover, the isoscaling phenomenon 
is actually observed experimentally for cold fragments. 
The sequential decay of hot primary fragments may not 
affect much the isoscaling coefficient. We notice that this 
is still a matter of hot debate depending on the model 
calculations 39, 40]. Therefore, with this view the small 
value of Csym for cold fragments extracted in the isoscal- 
ing experiments would indeed indicate that the fragments 
have dilute internal density. This would require a deep 
reconsideration of the statistical models from which Eq. 
© was derived. 

To this end, it is necessary to repeat here some remarks 
made by one of us at the WCI3 meeting^Jl- While the 
Eq.® is a good approximation within the grand canon- 
ical statistical model for multifragmentation, the isoscal- 
ing coefficient a is sensitive to the density dependence 
of the symmetry energy not because of the Csym which 
should be the symmetry energy of normal nuclear matter, 
but rather because of the Z/A ratios of the fragmentating 
sources through the dynamical isospin fractionationl42j 
in the early stage of the reaction. This was also pointed 
out recently in ref.'i^. Unfortunately, since the Z/A ra- 
tio of the effective fragmentating source, if exists at all, 
is not directly accessible experimentally, simplified as- 
sumptions are normally made in the data analyses within 
both statistical and dynamical models. The efforts are 
normally misplaced on extracting the Csym as if it is the 
mentation, some experts take the view tll0,t tllG Csym 

is one depending on the density and temperature. Conse- 



temperature at sub-saturation densities. It is seen that 
in the temperature range considered here, the symme- 
try energy does not change much with temperature at 
a given density. It is interesting to mention that using 
the same Gogny interactions corresponding to x = and 
X — —1 [s^ l in the finite temperature HE approach, all 
of the major results discussed above are qualitatively re- 
produced '28]. The major physical features given by the 
MLZ thermal model are thus rather general. 



III. ISOSCALING IN HEAVY-ION COLLISIONS 

It has been observed in many types of reactions that 
the ratio R2i{N,Z) of yields of a fragment with proton 
number Z and neutron number N from two reactions 
reaching about the same temperature T satisfies an ex- 
ponential relationship R2i{N,Z) cx exp(a7V) 0, 0, 0, 
iSmSllilsiHiEl- In several statistical and dy- 
namical models under some assumptions [s^ |3 H^i it 
has been shown that the scaling coefficient a is related 
to the symmetry energy Csym {p, T) via 



■^CsymjP-, T) 

T 



A [{Z/Af 



(6) 



where ^[{Z/Af] = {Z^j A^f -{Z^l A^f is the difference 
between the {Zj A)"^ values of the two fragmenting sources 
created in the two reactions. 

Before we proceed, a few comments and discussions re- 
garding the validity of Eq. © and the physical meaning 
of Csym are in order. First of all, we notice that Eq. © 
is an approximation in equilibrium models and an em- 
pirical assumption in dynamical models where isoscaling 
is observed in generated events. Because of the different 
assumptions used in the various derivations, the validity 
of this equation is still disputable as to whether and when 
the Csym is actually the symmetry energy or the symme- 
try free energy. Moreover, the physical interpretation of 
the Csym{p,T) is also not clear, sometimes even contra- 
dictory, in the literature. The main issue is whether the 
Csym measures the symmetry energy of the fragmentat- 
ing source or that of the fragments formed at freeze-out. 
This ambiguity is also due to the fact that the derivation 
of Eq. JBJ is not unique. In particular, within the grand 
canonical statistical model for multifragmentation |33ll3^ 
the Csym refers to the symmetry energy of primary frag- 
ments. While within the sequential Weisskopf model in 
the grand canonical limit |33j| it refers to the symmetry 
energy of the emission source. 

Based on the grand canonical model for multifrag- 
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quently, the extracted Cgym may thus contain some in- 
formation about the density dependence of the symmetry 
energy that is actually carried by the Z/A ratios of the 
fragmentating sources. 

On the other hand, within the sequential Weisskopf 
model in the grand canonical limit |lq . all quantities in 
Eq. refer to the emission source. In particular, the 
Caym itsclf in Eq. © reflects the bulk symmetry en- 
ergy of the low density fragmentating source. Besides, 
because of the way by which the data is analyzed as we 
mentioned above, the experimentally extracted Cgym also 
contains information about the symmetry energy through 
the Z/A ratios of the fragmentatig sources. Here we thus 
use broadly a working assumption that the Csym reflects 
the symmetry energy of bulk nuclear matter. We notice 
that this assumption can be justified only in the sequen- 
tial weisskopf model. Nevertheless, it is interesting and 
reassuring to note that this assumption is consistent with 
the statement in Ref. Ise| that the Csym is the symmetry 
energy of uniform nuclear matter at a reduced density. 
Within this picture it is natural for the Csym to have 
values smaller than the symmetry energy of normal nu- 
clear matter. One would thus have no difficulty with the 
Fisher hypothesis for fragment formation. At freeze-out, 
fragments are formed at their normal density in a large 
volume. Only the average density in the freeze-out vol- 
ume is small. 

We would also like to stress here that the Cgym {p, T) 
extracted from studying the isoscaling coefficient a is the 
total symmetry energy at the finite temperature T. It 
should be distinguished from the symmetry energy at 
zero temperature extracted from transport model anal- 
yses of dynamical observables, such as the isospin diffu- 
sion, or that from studying the neutron-skins of heavy 
nuclei. In transport models, the nucleon mean field or 
effective interaction for cold nuclear matter is used as an 
input. The zero temperature symmetry energy as a func- 
tion of density can thus be constructed analytically from 
the particular mean fields or effective interactions used 
in the calculations. Special cares should thus be taken 
when comparing the density functions for nuclear sym- 
metry energy extracted from different approaches and/or 
results from the same approach but for reactions at dif- 
ferent temperatures, especially at low densities. 

Within our working assumption that the extracted 
Csym from isoscaling analyses is the symmetry energy 
of the fragmenting source, we can then compare the ex- 
perimental Csym with our calculations of the symme- 
try energy for uniform nuclear matter. In Fig. 2, the 
experimental data taken from the SJY group at Texas 
A&M University (filled circles) [11 Ell and the INDRA- 
ALADIN collaboration at GSI (open squares) |20|,|23 are 
compared with the calculations. In the TAMU experi- 
ments, combinations of projectile-like fragments from pe- 
ripheral to semiperipheral collisions of 25 MeV/nucleon 
®^Kr and ^"^Ni beams on several neutron-rich targets were 
used in the isoscaling analyses. The symmetry energy 
was found to decrease quickly from about 25 MeV to 19 
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FIG. 2: (Color online) The relative evolution of symmetry 
energy as functions of temperature and density with two dif- 
ferent interaction-parts of the symmetry energy (see text). 



MeV as the temperature increases from about 4.8 MeV 
to 5.8 MeV HI. The INDRA@GSI data were obtained 
from the fragmentation of target-like residues following 
collisions of ^^C on -'^^^'^^^Sn targets at a beam energy of 
300 MeV/nucleon. The INDRA@GSI data indicate that 
the symmetry energy decreases from about 26 MeV to 16 
MeV as the temperature increases from about 6 MeV to 
9 MeV when the reaction goes from peripheral to central 
collisions poj . 

It is very interesting to compare the calculations us- 
ing both the X = (upper window) and the x = — 1 
(lower window) interactions with the experimental data. 
The comparison then allows us to estimate the required 
density of the fragment emitting source. At the respec- 
tively low temperatures reached in the peripheral reac- 
tions at TAMU and GSI, the calculated results indicate 
that the fragments are emitted from sources at densities 
only slightly below pQ. In the peripheral reactions, either 
the projectile-like or target-like residue is only slightly 
excited with little expansion. While at the higher tem- 
peratures reached in the more central reactions, the frag- 
ments are emitted from significantly diluted sources with 
densities depending on the interaction used. This pic- 
ture is consistent with dynamical model calculations of 
nuclear multifragmentations in the energy range consid- 
ered. More interestingly, it is seen that the experimen- 
tally observed evolution of the symmetry energy is mainly 
due to the change in density rather than temperature. 
Around the typical freeze-out temperatures reached in 
both the TAMU and the GSI experiments, the evolution 
of the symmetry energy due to the change in tempera- 
ture at a given density is rather small. It implies that the 
underlying origin of the observed decrease of the symme- 
try energy with the apparently increasing temperature 
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or centrality in these experiments is actually due to the 
accompanying decrease in freeze-out density. Comparing 
the two sets of data, it is seen that they actually paral- 
lel with each other in the common density range. Since 
the evolution of the symmetry energy is essentially inde- 
pendent of the temperature for the experiments consid- 
ered, the two sets of date thus indicate the same density- 
dependence of the symmetry energy as one expects. On 
the other hand, within the view that the extracted Cgym 
reflects the symmetry energy of the fragments at freeze- 
out, the two data sets are incompatible. 

From the above discussions, we can see that the evo- 
lution of the symmetry energy can be useful in explor- 
ing the density dependence of the interaction part of the 
nuclear symmetry energy. The latter is most uncertain 
but very important for many interesting questions in as- 
trophysics. With the interaction labeled a; = 0, the 
hottest point requires an average freeze-out density of 
about 0.62po and 0.49po for the TAMU-SJY and the 
INDRA@GSI data, respectively. While using the in- 
teraction labeled x = —1, the corresponding average 
freeze-out density is about O.Spo and O.GSpo, respec- 
tively. Therefore, the required freeze-out density de- 
pends strongly on the interaction part of the symme- 
try energy. The effect of using interactions of x = 
and .T -1 is about 30% and 40% for the TAMU and 
the INDRAQGSI experiment, respectively. Of course, 
for the purpose of extracting the interaction part of the 
symmetry energy from the isoscaling analyses, it is nec- 
essary to know not only the freeze-out temperature but 
also the density when the fragments are emitted from 
the reactions. Fortunately, inferring both the freeze-out 
temperature and density in isoscaling analyses has been 



shown feasible very recently by the TAMU-JBN group 
|3ll |. albeit largely based on model calculations. Since 
an independent determination of the freeze-out density 
will be very useful, observables known to be sensitive 
to the freeze-out density, such as the fragment correla- 
tion functions l45l | or source functions from imaging 
technig ues j46l 14^ . may be explored together with the 
isoscaling analyses. 



IV. SUMMARY 

In summary, within the degenerate Fermi gas model 
of Mekjian, Lee and Zamick it is shown that the experi- 
mentally observed evolution of the symmetry energy can 
be well understood. Furthermore, the evolution is found 
to be mainly due to the variation in the freeze-out den- 
sity rather than temperature in the reactions carried out 
under different conditions. The isoscaling analyses are 
thus useful for probing the interaction part of the nu- 
clear symmetry energy, provided that both the freeze-out 
temperature and density can also be inferred. 
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